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THE POTENTIAL FUNCTION AND 
FIELD OF FORCE OF A SPHEROID 
by 

James P. Murphy 
SUMMARY 

General formulas are presented for the potential function of a 
spheroid in rectangular coordinates. The components of the force, 
the partial derivatives of the potential with respect to the coordi- 
nates, are also derived. 
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THE POTENTIAL FUNCTION AND 
FIELD OF FORCE OF A SPHEROID 


INTRODUCTION 

With progress in tracking systems and the availability of more and better 
tracking data from the many artificial satellites launched and about to be launched, 
the need arises to include smaller forces in analyzing the motion of the satellites. 
Among these smaller forces are the so called "higher harmonics" in the earth's 
gravitational potential. The purpose of this report is to present general expres- 
sions depending only onm, the order, andn, the degree, for the components of 
the force due to a particular harmonic. With such a formulation in a system for 
analyzing data, generating tracking predictions, deriving harmonic coefficients 
etc., there would be no need to revise the programs to include any additional 
harmonic, but simply to imput the various (m, n) that one wished to consider. 


THE POTENTIAL FUNCTION 

The potential of a spheroid at a point, p, with spherical coordinates (r, 0, k) 
is 


co n 

» ■ tWLLW P n ( m ) (sin0) (C n ^ m ^ cosmA. + S n ^ m ^ sinmX) 


(1) 


n=2 m=0 


where /x is the product of the gravitational constant, G, times the mass of the 
spheroid, M, and b is the mean radius of the spheroid. The Legendre associated 
function, P n (m) (sin 0), is defined by 


P n (m) (t) 



(1 _ t 2 ) m /2 


d n+m ( t 2 _ 

dt n+m 


(2) 


The quantity R = U - n/r , is called the disturbing function and R n (m) is that part 
of the disturbing function due to the harmonic of order m and degree n. Therefore, 

t 

R ~ P ( m ) ( sin 0) ( C cos mk + S sin mk ) . (3) 

n r n+ln v ^ / \n n / 


k\ 


1 



One of the factors of R n (m) , P n (m) (sin 4>), will now be simplified: 


p n < m > (sinc£) - y (l-sin 2 0) 


m/2 


d n+m (sin 2 4>~ l) n 


d(sin 4>) 


n+m 


2 n n! 


cos m <fr 


d n+m (-1) 1 (") sin 2 < n_1 > 0j 


d( sin cp) 


n+m 


2 n n! 


(n-m)/2 [(n+m) even] 

(n-m-l_)/2 [(n+m) odd] r 2 ( n - i )1 1 

cos m 4> / C" 1 ) 1 (j) (n-m- 2i)! sin™' 2 * ' 

i = 0 


(4) 


After substituting Equation 4 into Equation 3, and after making use of the follow- 
ing relationships between spherical and rectangular coordinates: 


z 

sin <p - — 


m/2 (m even) 

( m — 1 )/ 2(m odd) 

cos m 0cosmA. - r“ m y 1 (~l) j ^2j) xm ” 2 ^ y2j 

j“0 


(m-2)/2 (m even) 

( m — 1 )/2 (m odd ) 

cos” <£ sin mA. = r“ m ^ (~l) k (2^+ x) x m ~ 2k-1 y 2k+1 

k=0 


where 

r = -^x 2 + y 2 + z 2 , 
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Equation 3 becomes 


R < m > 


Aib n 
2 n n! 


r (n+m+1) y(m, n) [C n < ra > a(m) + S n < m ) /3(m)] 


(6) 


where 


m/2 (m even ) 

(m- 1 )/2 (m odd) 

<*0) = r ( -iy ( 2 m j) if m/0 

j=0 


(m-2)/2 (m even) 

(m-l)/2 (m odd) 

/3(m) = 2] (-l) k ( 2k m + i) x m_2k_1 y 2 ^ if m / 0 

k=0 


7(m, n) 


(n-m)/2 (m + n even) 
( n — m - 1 ) / 2 (m+n odd) 

- £ 


i=0 


[2(n-i)]l 
(n - m - 2i ) ! 



n-m-2i 


(7) 


and 


a(0) = 1 if m = 0 

>8(0) = 0 if m - 0 . 


THE FIELD OF FORCE 

The field of force is given by the gradient of the potential function. Thus, 


F 


VU 


<?u r au r au- 

(?x i + dy i + dz ^ ’ 


( 8 ) 
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where I, ]*, and k are the orthogonal unit vectors in the spheroid centered co- 
ordinate system. If we denote any one of the coordinates x, y, or z by x p , then 



m, n 


(9) 


where 


<9R 

n 


r - (n+ m+l) 

2 n n! 



(n + m + 1) y(m, n) [C n (m) a(m) + S n < m) /3(m)] 


+ 


dy(m, n) 
? X p 


[C n < m > a(m) + S n ( m ) y 0(m)] 


+ y(m, n) 


C (m) 

n 


<?a(m) 

ax P 


+ s < m > 

n 


#) 

<9x 


( 10 ) 


The various partial derivatives of a, /3 , and y for x p equal to x , y , and z are now 
given. If m = 0, <9a(0)/<9x p = d/l(0)/<?x p = 0. If m f 0 


da(m) 

dx 


m/2 (m even) 

(m-l)/2 (m odd ) 

= ( -1 ) j (2j) (m- 2j ) x m_2j_1 y 2j 

TS 1 


<?a(m) 

dy 


m/2 (m even ) 

(m-l)/2 (m odd) 

2 L ] (_1)j (5) xm_2i y2i_1 

j=0 


da(m) 
d z ' 


0 


( 11 ) 
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*£(">) 

5x 


(m-2)/2 (m even)' 

(m — 1 )/2 (m odd) 

= ^ ( ~ 1)k (2* + 1) ( m ~ 2k ~ i)x m_2k ~ 2 y 2k+1 

k=0 


djB(m) 

dy 


(m-2 )/2 (m even ) 

(m- 1 )/2 (m odd) 

= L(-) k ( 2k m + i)( 2 k -)- 

k=0 


i-2k“l y 2k 


<^(m) _ ft 

~ - 0 • 


( 12 ) 


If m - n, dy(m, n)/dx p - 0. If m / n 


dy(m, n) 

<?X 


r(m, n) 


dy(m, n) 

?y 


r(m, n) 


dy(m, n) 

75 


x 2 + y 2 


(n-m-2)/2 (n+m even) 
(n-m-l)/2 (n+m odd) 


+ y 2 V~ ' . /n\ [ 2 ( n * )] ' /z\ 

r 3 / ( _1 )‘ \i) (n - m- 2i“ :r Ty! \ r J 


n-m“2i~l 


(13) 


where 


r(m, n) 


( n — m— 2 )/2 (n+m even) 


( n — m — 1 ) /2 (n+m odd) 

■ L 


i=0 


[2(n- i)]! 

(n - m - 2i ~ 1) ! 



n~m“2i 


(14) 


Therefore, by employing the relations given by Equations 8 through 14, the po- 
tential and field of force for any spheroid comprised of any combination of zonal 
(m = 0), tesseral (m / 0, n / m) or sectorial harmonics (m J 0,m = n) may be 
derived. 
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If the spheroid is rotating about the z-axis, then the influence of the tesseral 
and sectorial harmonics vary at a fixed point, p, in space due to the rotation. 

The force at p in a fixed, or inertial, coordinate system is given by 


_ (d u <9U \-» /<9U <JU \-, «?U - 

F - R 3 (-^)F - ^ cos a ~d^ sinaji + ^ sma +^cos ajj +J^k ,(15) 


where i ' , j ' , and k ' are the orthogonal unit vectors in the inertial system and 
a is the angular separation between T and X' . 


APPLICATION TO AN EARTH SATELLITE 


In the analysis of the motion of an earth satellite, many find it more con- 
venient to use the normalized harmonic coefficients C n (m) and S n (m) . They may 
be obtained from C (m) and S by 

n n ^ 


c (m) 

= N (m) C 

n 

n n 

s < m) 

n 

= n S 

n n 


(16) 


where 


N ( m ) 


N(2n +1) (n-m)! 
(n + m) ! 


- 1/2 


f N - 2 m / 0 

In = 1 m - 0 


(17) 


Values for many of the C n (m) and S n (m) for the earth may be found in Reference 1. 

For any earth satellite, the previously referred to quantity a becomes the 
Greenwich Mean Sidereal Time. Note that the rotation given by Equation 15 may 
be bypassed if only zonal harmonics are considered since the z coordinate is the 
same in both the rotating and inertial systems and r is invariant with respect to 
the rotation. 


Finally, for an axially symmetric potential (zonal harmonics only), the recom- 
mended potential function Reference 2 is 


U 


tL 

r 


1 " 2_]j n (7) P n ( Sint £) ' 

n= 1 _ 


(18) 
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Therefore, one must note that, if certain J n zonal harmonics are known but the 
tesseral and sectorial harmonic are to be considered also, in the equations above 
we have 


C n,0 = - J 

S n,0 = 0 . 
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